The chromatin landscape of the euryarchaeon Haloferax volcanii

Background Archaea, together with Bacteria, represent the two main divisions of life on Earth, with many of the defining characteristics of the more complex eukaryotes tracing their origin to evolutionary innovations first made in their archaeal ancestors. One of the most notable such features is nucleosomal chromatin, although archaeal histones and chromatin differ significantly from those of eukaryotes, not all archaea possess histones and it is not clear if histones are a main packaging component for all that do. Despite increased interest in archaeal chromatin in recent years, its properties have been little studied using genomic tools. Results Here, we adapt the ATAC-seq assay to archaea and use it to map the accessible landscape of the genome of the euryarchaeote Haloferax volcanii. We integrate the resulting datasets with genome-wide maps of active transcription and single-stranded DNA (ssDNA) and find that while H. volcanii promoters exist in a preferentially accessible state, unlike most eukaryotes, modulation of transcriptional activity is not associated with changes in promoter accessibility. Applying orthogonal single-molecule footprinting methods, we quantify the absolute levels of physical protection of H. volcanii and find that Haloferax chromatin is similarly or only slightly more accessible, in aggregate, than that of eukaryotes. We also evaluate the degree of coordination of transcription within archaeal operons and make the unexpected observation that some CRISPR arrays are associated with highly prevalent ssDNA structures. Conclusions Our results provide the first comprehensive maps of chromatin accessibility and active transcription in Haloferax across conditions and thus a foundation for future functional studies of archaeal chromatin. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-023-03095-5.


Introduction
Life on earth is now understood to be divided into two deep fundamental clades -Archaea and Bacteria.Archaea were only discovered as a separate branch of the tree of life in the 1970s [1], yet it was noticed very early on that they share a number of common features with the more organizationally complex eukaryotes, especially in the organization of their information processing cellular machinery.Based on these similarities it was suggested that eukaryotes evolved from archaea [2,3], a view strengthened in the phylogenomic era [4], and eventually solidified with the discovery of archaeal lineages such as the Lokiarchaeota [5].Thus, we now know that many of the complex cellular features that characterize eukaryotes trace their origins to their archaeal ancestry [6,7].
One of the most notable such features is nucleosomal chromatin.Nearly all eukaryotic genomes are packaged by nucleosomes, consisting of two tetramers of the four core histones H2A, H2B, H3, and H4, wrapping around ∼147 bp of DNA.These proteins are, with very rare exceptions [8,9], the most evolutionarily conserved among eukaryotes [10], in large part because aside from their packaging function they are also subject to a large number of precisely regulated posttranslational modifications (PTMs) at key residues [11], through which they play a pivotal role in all aspects of chromatin biology (transcription and its regulation, DNA replication, DNA repairs, mitosis, and others).
As early as the 1980s, it was noticed that some archaea possess proteins and structures similar to eukaryotic histones and nucleosomes [12,13].We now know that most archaea have histones genes [14][15][16] and that these histones are ancestral to the eukaryotic histones.Archaeal histones differ substantially from those in eukaryotes -while they share the core histone fold domain, they usually do not have the unstructured tails of H2A/H2B/H3/H4 that are the main sites of key PTMs.Archaeal histones also do not form octameric nucleosomes; instead, only one or a very small number of histone genes are found in archaeal genomes, and the structures they form are very different from those of eukaryotes.The diversity of histone sequences across the whole archaeal phylogeny is very large and still largely unexplored experimentally, but the available structural [17], biochemical, and modeling work suggests that in at least some species histones can form so-called hypernucleosomes or archaeasomes, consisting of a protein core of individual histones stacked next to each other, around which DNA is wrapped [15,18], ranging from 60 to 500 bp [19].It has also been proposed that archaeal histones exhibit an inherently dynamic association with DNA [19], in contrast to the much more stable association of nucleosomes with DNA in eukaryotes.
However, chromatin is not nucleosomal in all archaeal lineages -some lack histone genes and use other proteins such as Alba/Sac10b, Sul7d, Cren7, and CC1 [20,21] to package their genomes, and even for the ones that do have histone genes it is not always clear that histone proteins are in fact the main packaging component of chromatin (see the "Discussion" section further below).The latter situation is particularly notable in haloarchaea [22][23][24][25][26][27] -the tentative conclusion that their histones do not play a dominant role in DNA packaging is based on the low abundance of the HstA protein in Haloferax proteomics datasets [26,27] and on data suggesting that its functions are more analogous to those of a transcription factor [24].
Despite the relevance of archaeal chromatin to understanding the deep evolution of chromatin organization, up until now, the structure of archaeal chromatin has received little direct experimental investigation using modern genomic tools, with the exception of early MNase-seq studies nearly a decade ago that mapped nucleosomal positioning in the euryarchaeotes Haloferax volcanii [28] and Methanothermobacter thermautotrophicus and Thermococcus kodakarensis [29], and more recent MNAse-seq studies in Methanothermus fervidus [30], Thermoplasma acidophilum [31], and Methanocaldococcus jannaschii [32].Furthermore, very little is known about the relationship between chromatin structure and the regulation of gene expression in these organisms.
In order to begin to fill these gaps in our understanding of the organization of archaeal chromatin, we mapped chromatin accessibility and active transcription in Haloferax volcanii using a combination of bulk and single-molecule techniques such as ATACseq [33] (Assay for Transposase-Accessible Chromatin using sequencing), NOMe-seq/ dSMF [34] (Nucleosome Occupancy and Methylome sequencing/dual Single-Molecule footprinting) and KAS-seq [35] (kethoxal-assisted single-stranded DNA sequencing).We find that chromatin in H. volcanii exhibits similar features to that of eukaryotes on a broad level, with preferentially accessible promoter regions.However, unlike in eukaryotes, chromatin accessibility at promoters does not relate to transcriptional activity.Using single-molecule footprinting we estimate absolute protein occupancy levels over the H. volcanii TSSs to be comparable to, or possibly slightly lower than those in eukaryotes.However, unlike what is seen in most eukaryotes, we do not observe stably positioned nucleosome protection footprints, but rather only statistically elevated accessibility around promoters.We also revisit the question about the degree of coordination of transcriptional activity and chromatin accessibility within Haloferax operons and make the unexpected discovery that some CRISPR arrays are associated with very strong ssDNA signatures.

ATAC-seq reveals the open chromatin landscape of H. volcanii
In order to study chromatin accessibility in archaea, we adapted the ATAC-seq assay [33] to the Haloferax volcanii archaeon.H. volcanii is a halophile with a strong preference for very high salt concentrations in the growth medium (see the "Methods" section), which grows optimally at 42°C [36,37], and is a widely used archaeal model system.
The principle behind ATAC-seq is the very strong preference of the Tn5 transposase [33] for insertion into accessible DNA as opposed to tagmentation of protected (by nucleosomes, transcription factors, or other proteins) DNA.Tn5 insertion then tags accessible sites with landing sites for PCR primers, allowing for highly efficient amplification of open chromatin regions in the genome.
After extensive testing of a variety of different experimental protocols (fixation conditions and input cell numbers), we arrived at the following modifications of the standard ATAC protocol.First, because archaea are not eukaryotes and do not have a nucleus, we omitted the cell lysis and nuclei isolation step that is a standard feature of eukaryotic ATAC-seq protocols, such as the now standard omniATAC [38].Second, and most importantly, we reasoned that if the previously reported dynamic association with DNA of archaeal nucleosomes (or of other proteins) occurs in H. volcanii, optimal results might be obtained by introducing a crosslinking step into the standard ATAC protocol, which would "freeze" any histones and other structural proteins in place and not allow transposition into DNA that might change from protected to accessible during the duration of the transposition reaction.Indeed, comparing the TSS (transcription start site) enrichment generated without fixation and with light (0.1% formaldehyde) and strong (1% formaldehyde) fixation showed that stronger fixation produces higher TSS enrichment (Fig. 1A).We then compared the H. volcanii ATAC-seq TSS metaprofile with that from the previously published MNase-seq dataset and observed the expected inverse relationship (Fig. 1B).H. volcanii TSSs exhibit elevated accessibility in the 0 to − 400 bp upstream region, decreasing away from the TSS.
The H. volcanii genome consists of multiple replicons [39], with a main chromosome ("chr") and four plasmids of very different sizes -pHV4, pHV3, pHV1, and pHV2 (in order of decreasing size), which together comprise ∼30% of the total genome.To properly interpret sequencing data (which is typically normalized to total read coverage), we determined the relative copy number distribution of these replicons using tagmented naked genomic DNA control samples (Fig. 1C).The smallest plasmid -pHV2 -appears to exist in ∼26 copies for each main chromosome, pHV1 is found at ∼1.4 copies for each main chromosome, while the two large plasmids -pHV4 and pHV3exist in a 1:1 ratio to the main chromosome.
The H. volcanii ATAC-seq fragment length distribution is unimodal, peaking at 90-100 bp, and does not show an analog to the eukaryotic mono-, di-and tri-nucleosomal signature (Fig. 1D), even scaled down to the smaller protection footprint of archaeal histones.This suggests that Haloferax chromatin may not be packaged into abundant nucleosmal structures consisting of strings of closely positioned individual nucleosomes.
Peak calling using MACS2 [40] (Fig. 1E) and genome browser visualization of ATACseq profiles (Fig. 1F-G) revealed an accessibility landscape largely reminiscent of that in eukaryotes with compact genomes such as the budding yeast Saccharomyces cerevisiae [41] -nearly all ATAC-seq peaks are located within 200 bp of an annotated TSS, and these peaks are very strong and localized.
As previous studies of chromatin openness in bacteria have reported the existence of very large domains of lower and higher accessibility [42], we wondered whether the same is observed in Haloferax.We do not observe such domains in our datasets (Fig. 1I).Recently, an ATAC-seq dataset was reported from the crenarchaeote Sulfolobus islandicus [43], which lacks histones, and instead packages its genome mainly through Alba/ Sac10b proteins [44,45].In that species, large domains similar to those in bacteria were reported.This observation may suggest that such large-scale domains of elevated chromatin accessibility are a feature associated with the lack of nucleosomal chromatin in prokaryotes, while archaea that contain histone genes, such as H. volcanii, exhibit eukaryote-like organization, but such a generalization is contingent on Haloferax chromatin being in fact nucleosomal.We also reexamined the Sulfolobus islandicus dataset and found that it displays a much more modest TSS enrichment than that seen in H. volcanii, which is also more narrowly concentrated around the TSS position (Additional file 1: Fig. S1).
Finally, we observed an anti-correlation between ATAC-seq signal and genomic GC content (Additional file 1: Fig. S2).Haloferax volcanii exhibits a rather high average GC content of 65% [39], decreasing to 58% in intergenic regions, and areas with even higher GC content (> 65%) show markedly lower ATAC-seq signal.This observation is corroborated by the available external MNase-seq dataset, which shows a positive correlation with GC content (i.e., the inverse of ATAC-seq, as expected) and naked DNA controls (which show no correlation with GC content, indicating that PCR biases during sequencing library preparation are not the reason for the observed patterns).

Absolute DNA occupancy/protection levels in H. volcanii
While ATAC-seq is immensely helpful for identifying the location of accessible regions in the genomes and measuring their relative accessibility, it is a bulk method that does not provide information about the absolute levels of protection/accessibility in the genome.Instead, absolute accessibility must be measured by either restriction digestion-based or enzymatic labeling single-molecule methods.To quantify absolute occupancy/protection levels in the H. volcanii genome, we applied NOMe-seq [46] and dSMF [34] to Haloferax chromatin.These methods rely on the preferential methylation of accessible cytosine nucleotides (5mC) by a recombinant methyltransferase that modifies specifically in GpC contexts (NOMe-seq) or a combination of methyltransferases that label both GpC and CpG (dSMF).
However, these methods are potentially confounded by the presence of endogenous methylation in either context.Fortunately, in the case of H. volcanii endogenous DNA modifications have been previously studied using PacBio single molecule sequencing, and no CpG and GpC modifications were found.Instead, only two restriction modification system-associated modifications in different contexts, specifically, 4-methylcytosine in a C(m4)TAG context and N6-methyladenine in a GCA(m6)BN6VTGC context were identified [47].
Figure 2A-C show the metaprofiles of average methylation around H. volcanii TSSs for NOMe-seq and dSMF datasets generated from exponentially growing and stationary cultures (post log-phase in the growth curve).We observe baseline absolute protection levels around 84-85% in the exponentially growing cells and ∼89% in stationary cells.For comparison, analogous studies in eukaryotes, such as the budding yeast S. cerevisiae [41,48,49], have shown absolute protection levels around 90% (± 5%).Thus, Haloferax chromatin exhibits broadly similar, though perhaps somewhat lower levels of protection than what is observed in conventional eukaryotes.
The base-pair resolved nature of these single-molecule methods enabled an observation of a feature not readily apparent in ATAC-seq and MNase-seq datasets -a protection footprint immediately upstream of the TSS.We also observe this feature as a protection footprint in a few percent of single molecules at individual promoters (Fig. 2D).At present we are not able to confidently identify its functional association -its width is likely too small for it to be a positioned -1 nucleosome, and we hypothesize it may correspond to one of the complexes involved in the archaeal transcriptional cycle, analogous to how similar protection footprints associated with the RNA polymerase and the preinitiation complex (PIC) in eukaryotes have been observed in dSMF datasets [34].
On the other hand, unlike this unique protection footprint, we do not observe strongly positioned individual nucleosome-like features along the Haloferax genome, such as those seen in conventional eukaryotes (Fig. 2D).However, our NOMeseq and dSMF datasets were sequenced at an effective depth of ∼100 × for fragments of width 200 bp.To determine if higher sequencing depth would clearly reveal such putative positioned nucleosomes (or other proteins), we turned to the pHV2 plasmid, which, as previously discussed, exists in high copy numbers in H. volcanii cells.Over the pHV2 plasmid, we obtained ∼200 × coverage for fragments of length 250 bp (Fig. 2E) and ∼1200 × coverage for fragments of length 200 bp (Fig. 2F).These high-depth maps also reveal considerable heterogeneity of footprints and accessible sites.These observations are not consistent with the existence of abundant strongly positioned nucleosomes (or other structural proteins) associated with the Haloferax genome; they can be explained by dynamic and/ or positionally unstable association of chromatin proteins with DNA.
Finally, we also note that the absence of distinct states in the Haloferax NOMe-seq/ dSMF data points to all copies of its various replicons existing in mostly the same state.The two small Haloferax plasmids are certainly polyploid, and highly so in the case of pHV2, but quite likely even the main chromosome exists in multiple copies in individual cells.We see no evidence that some of these copies may exist in a different chromatin state (e.g., silent versus active) than the others.

The ssDNA and active transcription landscape in the H. volcanii genome
We then turned our attention to the landscape of active transcription in H. volcanii.To this end, we used KAS-seq [35], which measures with high specificity the presence of single-stranded DNA in the genome, by labeling unpaired guanine nucleotides with N 3 -kethoxal, using click-chemistry to add a biotin moiety, then enriching these fragments via a streptavidin pulldown.Most ssDNA is usually found within the transcriptional bubbles associated with RNA polymerase molecules engaged with DNA [35].KAS-seq provides several advantages in the H. volcanii context.First, due to the absence of straightforward methods to deplete H. volcanii ribosomal RNA (rRNA) from RNA sequencing libraries analogous to polyA-selection in eukaryotes [50,51], it enables measurement of transcriptional activity at a much lower cost than deep RNAseq experiments.Second, and most importantly, it measures actively engaged polymerase molecules, and thus unlike the steady-state transcript levels that conventional RNA-seq quantifies, it provides a means to quantify active transcriptional activity.Third, it also identifies other ssDNA structures, such as those resulting from paused polymerase molecules, G-quadruplexes, and others.
We carried out a time course analysis of Haloferax growth and applied both KASseq and ATAC-seq during the "exponential" log-phase of growth, and on the "stationary" post-log phase stage, as well as on "standing" cultures, which had been left at room temperature for ∼1 week.We also carried out KAS-seq on exponentially growing cells that were then incubated at different temperatures -the typical growing temperature of 42 °C, 37 °C, 23 °C and a cold shock at 4 °C for 4 h.
At a global level, we observe uniform levels of KAS signal along the length of H. volcanii chromosomes (Fig. 3A), with sharp localized peaks.KAS-seq measurements in H. volcanii are highly reproducible between experimental replicates (Fig. 3B).Locally, at the level of individual genes, we observe a combination of high-signal peaks at the promoters of some genes and elevated KAS-seq signal along gene bodies (Fig. 3C-D).This feature is also in KAS-seq metaprofiles over all genes (Fig. 3E), suggesting that in H. volcanii RNA polymerases spend substantial amount of time associated with the TSS, perhaps in a paused state analogous to that observed in metazoans [52] or the archaeon Sulfolobus solfataricus [53].

CRISPR arrays
In the process of optimization of the KAS-seq assay in Haloferax, we carried out KAS-seq on a H. volcanii culture that had been left standing at room temperature for ∼3 months.These data revealed that CRISPR arrays can become highly enriched for KAS-seq signals in these viable, yet dormant, cultures.CRISPR (clustered regularly interspaced short palindromic repeats) arrays are a key element in the defense systems against foreign genetic material of many prokaryotes and consist of multiple identical repeats interspersed with non-repetitive sequences that target foreign plasmids and phages, together with a set of Cas genes.H. volcanii is one of the prokaryotic systems where these elements were first originally observed [54][55][56].
In standing Haloferax cultures, transcriptional activity is likely suppressed, as the cells enter a dormant state.Consistent with a transcriptionally inactive state, we observe largely flat, low levels of KAS-seq signal from this standing culture KAS-seq dataset (Fig. 4A).However, we also observed a single large, sharp peak of KAS-seq signal on the pHV4 plasmid.This peak resides between the second (as numbered in the available genome annotation) CRISPR array in the H. volcanii genome [57] and its associated Cas6 gene (Fig. 4B); a previously reported small RNA (sRNA) -s479 [58] -is also located in between Cas6 and the CRISPR array, but our observed KAS-seq peak is not associated with this putative promoter, but is instead situated downstream of the array.While this KAS-seq signal peak is also found in all other conditions we assessed, it stands out in the long-term standing culture due to the absence of other peaks that result from the active transcription of other genes.
Curiously, only the second CRISPR array in H. volcanii displays this strong ssDNA structure, while the other two do not (this is true across all assayed conditions).However, all three arrays show elevated chromatin accessibility in ATAC-seq datasets -an accessibility signal that is not focused on the beginning of the array but covers its whole length (Fig. 4C).Possible interpretations of these observations are discussed below.

Coordination between chromatin accessibility and transcriptional activity within H. volcanii operons
Being prokaryotes, archaea often have genes organized into operons [59], with multiple genes transcribed as a single unit, which are therefore expected to share a common promoter and to exhibit similar levels of active transcription.However, transcription of these operons is still little studied using modern genomic tools.To address this gap, we used our KAS-seq and ATAC-seq data, which provide information about the chromatin accessibility and active transcription in different conditions, to investigate the extent of coordination between the transcriptional activity of different units in operons in the H. volcanii genome.
We first inspected the two H. volcanii operons that include rRNA genes [60]. Figure 5A shows KAS-seq and ATAC-seq profiles along rRNA genes in the different conditions we Fig. 4 Abundant ssDNA structures associated with some H. volcanii CRISPR arrays in specific conditions.A Global KAS-seq profiles over each of the five H. volcanii chromosomes in a long-standing culture (∼3 months) reveal an extremely strong ssDNA peak associated with one of the CRISPR arrays on the pHV4 plasmid.B KAS-seq signal levels around pHV4 plasmid CRISPR arrays in different conditions.C KAS-seq and ATAC-seq levels around all three H. volcanii CRISPR arrays in different conditions assayed.We observe a largely uniform KAS-seq profile in exponentially growing cells, and generally elevated chromatin accessibility (which might be associated with very active transcription); a more non-uniform pattern is seen in cold-shocked cells kept at 4 °C, where lower KAS-seq levels are seen over the large subunit (LSU) rRNA relative to the small subunit (SSU), with various intermediate states in other conditions.
More interesting patterns are seen in operons comprised of protein-coding genes.Figure 5B shows a gene array consisting of A-type ATP synthase subunits, for which distinct KAS-seq peaks are seen at the beginning of the operon as well as in between genes in the middle of the operon.Furthermore, KAS-seq levels are not uniform over the gene bodies of all genes.
We examined multiple other operons (Fig. 5C-D) and Additional file 1: Fig. S3, which reveal a very diverse picture of the extent of coordination between the transcriptional activity over individual genes within an operon -internal operon peaks are observed for multiple operons, while there are also other operons where KAS-seq signal is more uniform.
In some cases (e.g., Fig. 5C), these internal KAS-seq peaks are also associated with matched ATAC-seq peaks.Thus, one interpretation of these observations is that not all these operons are true operons (even though they consist of functionally related genes), but instead independent initiation and regulation of transcription from internal TSSs may be occurring.This interpretation is particularly supported in the cases where ATAC-seq peaks are seen at the beginning of genes, and where the baseline gene-body KAS-seq signal differs greatly between different sections of the operon and is also consistent with previous observations of internal promoters inside archaeal operons based on transcript 5′-end mapping [50,[61][62][63].On the other hand, internal KAS-seq peaks might also arise from very strong and immediate coupling between transcription and translation, i.e., if the process of initiation of translation at internal positions in the operon somehow leads to the polymerase pausing at certain sites.A, A, B, C, D, E, F, I, K, and H. C RNA polymerase II subunits.D NADH dehydrogenase-like complex subunits A, B, CD, H, I, J1, J2, K, L, M, and N

Chromatin accessibility does not correlate with transcriptional activity in H. volcanii
In eukaryotes, the regulation of chromatin accessibility at regulatory elements (promoters and enhancers) is key to gene regulation, as nucleosomal chromatin is generally refractive to occupancy by regulatory proteins and to active transcription [64], and while perfect correlation between accessibility levels at promoters and gene expression is rare, open chromatin states are generally associated with increased transcriptional activity.
In contrast, the relationship between chromatin accessibility and transcriptional activity in archaea has not been systematically studied as chromatin accessibility has not been mapped globally, across conditions, and in conjunction with global measurements of active transcription.
We first identified differentially accessible promoter regions between the different conditions we studied (Fig. 6A).In contrast to an a priori expectation that changes in gene expression would be associated with shifts in chromatin accessibility levels around promoter regions, we did not find strong changes between exponentially growing and stationary cells (Fig. 6A).We observed large apparent differences in ATAC-seq signal in each of those two conditions and standing cultures (Additional file 1: Fig. S4), but in those comparisons, the profiles are highly skewed towards increased higher accessibility in the actively growing and stationary cells instead of showing the typical more symmetric changes (i.e., accessibility increasing and decreasing in both directions) between two conditions.
We believe this pattern is due to the dormant state of standing cultures in which we do not observe as strong ATAC-seq peaks as are observed in actively growing cells.In stark contrast to the lack of accessibility changes between non-dormant conditions, we find a large number of genes that display strong differential KAS-seq signals over their gene bodies (Fig. 6B).Thus while we observe no major changes in chromatin accessibility across these two conditions, we do note large-scale changes in RNA polymerase occupancy (as inferred by KAS-seq).
We then quantified the degree of correlation between KAS and ATAC signals (Fig. 6C-E) and found no significant correlation.We also found no correlation between the level of changes in chromatin accessibility and changes in RNA polymerase association with DNA (as inferred by KAS-seq) between exponential and stationary cells (Fig. 6F).
These global observations are supported by the study of individual loci.Figure 7A shows ATAC-seq levels and KAS-seq levels over all genes across different conditions, clearly indicating that these two signals are decoupled from one another.Figure 7B and C depict other examples of genes for which transcriptional activity shifts between conditions yet ATAC-seq profiles are largely identical.
We thus conclude that based on the currently available data, the modulation of chromatin accessibility does not appear to be a major determinant/correlate of transcriptional activity in Haloferax archaea.

Discussion
In this study, we adapted and applied methods for global profiling of chromatin accessibility and ssDNA in the euryarchaeote Haloferax volcanii, revealing the chromatin architecture of this representative of the haloarchaea.We identified several convergent and divergent characteristics with respect to those of conventional eukaryote properties.
The H. volcanii genome displays an accessibility landscape very similar to that of eukaryotes with compact genomes such as budding yeast -accessibility peaks are almost exclusively found very close to promoters.Absolute accessibility/protection levels are similar, perhaps slightly lower than those in budding yeast, with a baseline protection level of 85-90%, although these numbers need to be interpreted with some caveats.In yeast and other eukaryotes, applying NOMe-seq/dSMF to fixed and native chromatin returns comparable absolute protection values (unpublished data), but it is possible that in Haloferax chromatin is nevertheless a bit more "open" than in eukaryotes, as the fixation step provides a "frozen" in time snapshot of protein occupancy on DNA while association of proteins with DNA can still be more dynamic than it is for eukaryotic nucleosomes.We do not observe the strongly positioned nucleosome-like features such as those typical in eukaryotes in Haloferax, but instead observe a heterogeneous picture of footprints, consistent with such a dynamic association of chromatin proteins with DNA.
These observations do make a certain amount of sense in the light of recent reports rejecting the nucleosomal packaging of haloarchaeal genomes [22][23][24] but are also puzzling on their own.Once again, the picture of chromatin accessibility revealed by ATACseq in Haloferax is nearly identical to that of conventional eukaryotes with nucleosomal chromatin and densely packed genomes.At an arbitrary locus, one would note very few general differences between a Haloferax ATAC-seq signal track and a yeast ATACseq signal track.Both exhibit strong localized peaks at promoters, and similar levels of absolute protection/occupancy. What protein can provide such high levels of physical protection is at present unknown.Experiments involving in vitro reconstitution of the association of archaeal histones and other putative chromatin proteins with DNA as well as their heterologous expression in yeast might shed light onto this question in the future.
We also note that, unlike what has been reported about bacteria and archaea without histones, the H. volcanii genome does not exhibit large-scale domains of diminished and elevated accessibility.
In contrast to the norm in eukaryotes, accessibility at Haloferax promoters does not correlate with transcriptional activity.This intriguing observation will require further functional dissection in future work.The idea that chromatin accessibility is not necessary for transcription in archaea is supported by previous observations that transcription by the archaeal RNA Polymerase is slowed, but not blocked by archaeal nucleosomes [65].However, the molecular determinants of the observed accessibility remain unclear.Nearly all promoters in H. volcanii show some level of accessibility (Fig. 7A), but their levels differ greatly between individual genes.How these differential states are specified, and whether they might in fact change in conditions that we have not assayed remains to be determined.MNase-seq studies Methanothermobacter thermautotrophicus and Thermococcus kodakarensis [29] have suggested that nucleosome positioning in those organisms is significantly influenced by DNA sequence, but no such strong association was reported for Haloferax volcanii [28].In the currently available datasets we observe anti-correlation between chromatinization and high genomic GC content, but whether this is the primary determinant of nucleosomal or other protein occupancy, and whether this correlation can account for the large differences in promoter accessibility combined with a general absence of such strong peaks elsewhere in the genome observed in Haloferax remains uncertain.Generating chromatin accessibility and nucleosome positioning data across other archaeal clades and also in multiple closely related species will allow us to generalize these observations and to train fully powered models that relate sequence to chromatin accessibility, potentially identifying such determinants.
We also made the observation that operons in Haloferax display non-uniform levels of single-stranded DNA signal consistent with transcriptional activity, and may in fact consist of multiple distinct transcriptional units.The phenomenon of independent transcription of operon genes has been suggested by some previous studies in Haloferax and Sulfolobus [66][67][68], but this is again an observation and a hypothesis that needs to be generalized to and tested in not only more archaeal species, but also in bacteria, where the application of KAS-seq to the study of transcriptional activity may also result in unanticipated findings.In Haloferax we were only able to examine several dozen of unambiguous operons (e.g., unidirectional arrays of functionally related genes).
Bacteria are also highly relevant to the other surprising observation we made -the strong ssDNA structure present at the second Haloferax CRISPR array, especially in dormant cells that are otherwise mostly transcriptionally silent, but not at the other two CRISPR arrays.The second CRISPR array is uniquely associated with the Cas6 gene.Cas6 is the endoribonuclease that generates guide RNAs [57,69], thus one possible explanation for the strong ssDNA peak between the CRISPR array and Cas6 is that it represents paused RNA polymerase at the Cas6 promoter.If this explanation is correct, the fact that Cas6 is the only gene in the genome with this property in dormant cells is remarkable and perhaps points to the importance of retaining the ability to process CRISPR transcripts even in a dormant cellular state.Alternatively, the ssDNA structure might be related to the transcription of the CRISPR array itself; however, such an explanation does not speak to the uniqueness of the KAS-seq signal at the second CRISPR array and the absence of the same strong KAS peaks at the other two CRISPR arrays.The second CRISPR array is also separated from both the Cas6 gene and the KAS-seq peak by the s479 sRNA, the reported function of which is involved in the regulation of zinc transport [58]; how this second array might fit into the overall picture is not clear.The functional significance of elevated chromatin accessibility over CRISPR arrays is also currently unknown.As prokaryotes exhibit an immense variety of CRISPR systems and number and organization of CRISPR arrays, mapping these properties in multiple other prokaryotes would be highly informative.

Haloferax volcanii cell culture
H. volcanii cells were obtained from the DSMZ German Collection of Microorganisms and Cell Cultures GmbH (Cat # 3757) and cultured in Halobacterium media [70,71], prepared as follows: 7.50 g casamino acids, 10.00 g yeast extract, 3.00 g sodium citrate, 2.00 g KCl, 20.00 g MgSO 4 × 7 H 2 O, 0.05 g FeSO 4 × 7 H 2 O, 0.20 mg MnSO 4 × H 2 O, and 250.00 g NaCl were mixed with distilled water in a total volume of 1 L. The media were then autoclaved and allowed to cool.H. volcanii was typically grown at 42 °C, except where otherwise indicated.Cultures were stored at room temperature when not actively growing.

Haloferax volcanii genome assembly and annotations
For all analyses, the genome assembly and annotation for the Haloferax volcanii DS2 strain, downloaded from the NCBI database, and also matching the haloVolc1 version on the UCSC Microbial Genome Browser [72] (http:// micro bes.ucsc.edu/), was used.The UCSC Microbial Genome Browser was used for visualization of genome browser tracks.

ATAC-seq experiments
Several variations of the ATAC-seq assays were tested.As H. volcanii is an archaeon, i.e., a prokaryote without a nucleus, and as it does not have a cell wall (as many other prokaryotes do), the nuclei isolation step typical for ATAC-seq protocols used in eukaryotes was omitted.
For crosslinked ATAC-seq, cells were fixed by adding 37% formaldehyde (Sigma) at a final concentration of either 0.1% or 1% and incubating for 15 min at room temperature.Formaldehyde was then quenched using 2.5 M glycine at a final concentration of 0.25 M. Cells were subsequently centrifuged at 10,000 g for 2 min, washed once in 1 × PBS, and centrifuged again at 10,000 g for 2 min.Transposition was carried out as above for 15 min.The reaction was stopped with the addition of 150 μL IP Elution Buffer (1% SDS, 0.1 M NaHCO 3 ) and 2 μL Proteinase K (Promega, Cat # MC5005), then incubated at 65 °C overnight to reverse crosslinks.DNA was isolated by adding an equal volume of 25:24:1 phenol to chloroform to isoamyl solution, vortexing and centrifuging for 3 min at 14,000 rpm, then purifying the top aqueous phase using the MinElute PCR Purification Kit, eluting in 10 μL EB buffer.Libraries were generated as described above.

KAS-seq experiments
KAS-seq experiments were carried out following the previously published protocol [35] with some modifications in the sequencing library generation part.
Briefly, a 500-mM N 3 -kethoxal solution was brought to 37 °C and then added to 2 mL of culture at a final concentration of 5 μM.Cells were then incubated for 5 min at 37 °C in a ThermoMixer at 1000 rpm.
Cells were then pelleted by centrifugation at 10,000 g for 1 min, resuspended in 200 μL 1 × PBS buffer, and DNA was immediately isolated using the Monarch Genomic DNA Purification Kit (NEB, Cat # T3010S), with the modification that elution was carried out with 50 μL 25 mM K 3 BO 3 solution (pH 7.0).
Libraries were built on beads using the NEBNext Ultra II DNA Library Prep kit (NEB, Cat # E7645L).Biotin pull-down was initiated by pipetting 20 μL Dynabeads MyOne Streptavidin T1 beads (ThermoFisher Scientific, Cat # 65306) into DNA lobind tubes.Beads were separated on magnet, resuspended in 200 μL of 1 × TWB buffer (Tween Washing Buffer; 5 mM Tris-HCl pH 7.5; 0.5 mM EDTA; 1 M NaCl; 0.05% Tween 20), then separated on magnet again and resuspended in 300 μL of 2 × BB (Binding Buffer; 10 mM Tris-HCl pH 7.5, 1 mM EDTA; 2 M NaCl).The DNA (130 μL) was added together with 170 μL 0.1 × TE buffer, and incubated at RT on a rotator for ≥ 15 min.Beads were separated on a magnet, resuspended in 200 μL of 1 × TWB, and incubated at 55 °C in a Termomixer for 2 min with shaking at 1000 rpm.Beads were again separated on magnet and the 200-μL 55 °C TWB wash step was repeated.Beads were separated on a magnet and resuspended in 50 μL 0.1 × TE.
End repair was carried out by adding 7 μL NEB End Repair Buffer and 3 μL NEB End Repair Enzyme, incubating at 20 °C for 30 min, then at 65 °C for 30 min.
End repair was followed by adaptor ligation by adding 2.5 μL NEB Adaptor, 1 μL NEB Ligation Enhancer, and 30 μL NEB Ligation Mix, incubating at 20 °C for 20 min, then adding 3 μL USER Enzyme and incubating at 37 °C for 15 min.Beads were separated on a magnet, resuspended in 200 μL of 1 × TWB, then incubated at 55 °C in a Thermomixer for 2 min with shaking at 1000 rpm.Subsequently, beads were separated on a magnet and resuspended in 100 μL of 0.1 × TE, separated on a magnet again, resuspended in 15 μL of 0.1 × TE Buffer, and transferred to PCR tubes.
Beads were then incubated at 98 °C for 10 min, and libraries were amplified by adding 5 μL of i5 primer, 5 μL of i7 primer and 25 μL of 2 × Q5 Hot Start Polymerase Mix, using

Fig. 1
Fig. 1 Archaeal ATAC-seq and the open chromatin landscape of H. volcanii.A, B Adaptation and optimization of the ATAC-seq assay to the archaeal context.A Distribution of TSS ratio scores (see the "Methods" section for details) for native, 0.1%-and 1%-formaldehyde ATAC-seq libraries.C Fragment length distribution in H. volcanii ATAC-seq datasets.D Estimated relative copy number of H. volcanii chromosomes.Genomic DNA was tagmented and amplified (n = 4) and normalized read coverage was estimated for each chromosome/ plasmid.The average ratios are shown.E Distribution of MACS2 ATAC-seq peaks relative to TSSs.F, G Representative browser snapshots of ATAC-seq profiles along the H. volcanii genome.F High reproducibility of H. volcanii chromatin accessibility measurements using ATAC-seq.Shown is the between-replicate correlation over TSSs in RPM (reads per million) units.I Global ATAC-seq profile over each of the five H. volcanii chromosomes.The number in brackets corresponds to the magnification of the true proportional size of plasmids relative to the main chromosome

Fig. 2
Fig. 2 Absolute DNA occupancy/protection levels in H. volcanii.A-C TSS metaprofiles in different conditions (two replicates of an exponentially dividing culture, and a stationary culture).D Single-molecule map (250-bp) around a main-chromosome TSS.Black indicates unmethylated and therefore protected sites, and gray indicates methylated and thus accessible sites.E, F Single-molecule maps over the pHV2 plasmid: 250-bp window map (E) and very high coverage (≥ 1200 single molecules) 200-bp window map (F)

Fig. 3
Fig. 3 The ssDNA and active transcription landscape in the H. volcanii genome as measured by KASseq.A Global KAS-seq profiles over each of the five H. volcanii chromosomes in an exponential culture.B High reproducibility of active transcription measurements using KAS-seq.C-D Representative browser snapshots of KAS-seq profiles along the H. volcanii genome.E KAS-seq metaprofile along H. volcanii gene bodies

Fig. 5
Fig. 5 Coordination between chromatin accessibility and transcriptional activity within H. volcanii operons.The black bar shows the operon boundaries.A Ribosomal RNA operon.Note that the tracks shown here were generated by including multimapping reads (see the "Methods" section for details).B A-type ATP synthase subunitsA, A, B, C, D, E, F, I, K, and H. C RNA polymerase II subunits.D NADH dehydrogenase-like complex subunits A, B, CD, H, I, J1, J2, K, L, M, and N

Fig. 6
Fig. 6 Chromatin accessibility does not correlate with transcriptional activity in H. volcanii.A Differential chromatin accessibility between exponential and stationary conditions.B Differential KAS-seq levels between exponential and stationary conditions.C-E Lack of correlation between KAS and ATAC signals in exponential and stationary conditions.F Lack of correlation between changes in chromatin accessibility and changes in transcriptional activity

Fig. 7
Fig. 7 Chromatin accessibility does not correlate with transcriptional activity in H. volcanii.A Genome-wide heatmaps of ATAC-seq and KAS-seq signals around H. volcanii TSSs, sorted by KAS-seq levels in the exponential condition.B, C Representative snapshots of genes with significantly altered transcriptional activity between the exponential and standing conditions, but no corresponding changes in chromatin accessibility